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g Abstract 

. The baryon-number density formed in relativistic nuclear collisions, versus the chemical potential 

Q-i' of the freeze-out states, is systematically studied on the basis of existing measurements. A remark- 

■ able power-law behaviour of the baryon-number susceptibility is found at the SPS, consistent with 

m 

the existence of a QCD critical point at ^b,c — 214 MeV, Tc ~ 155 MeV. The equation of state in 

'q^, different asymptotic regimes of the critical region is also examined and confronted with freeze-out 

Qh' states in these experiments. 
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Baryons produced in collisions of nuclei at relativistic energies play important role 
in the identification of critical states in strongly interacting matter The baryon- 

number density nsix) is an appropriate order parameter of the critical system, equivalent 
to the chiral field since they both incorporate in their behaviour, the universal power 
laws and the related exponents of the QCD critical point [1]. A complete set of observables 
associated either with the baryon-number density or the sigma field, may form together with 
the corresponding measurements in experiments with nuclei, the basic ingredients and tools 
for the construction of an observational theory of critical QCD matter [2|,0]- A search for 
power-law fluctuations in the isoscalar component of tt^tt" pairs near the two-pion threshold 

[in 

or in the net-baryon system produced at midrapidity, are examples of such observables |2j,p|- 
In this Letter we search for a singular behaviour of baryon-number susceptibility xb 
exploiting the existing measurements of baryochemical potential in a series of freeze-out 
states ranging from SPS to RHIC energies. The guidelines in this investigation come from 
the universal properties of the equation of state for 3d Ising systems which specify the 
universality class of the QCD critical point {4]. The observables in this treatment are (a) 
the baryon number density Ub, properly averaged in a space-time region occupied by the 
system and (b) the baryochemical potential fiB and the temperature T of the corresponding 
freeze-out state. The behaviour near the critical point is described by the order parameter 
Pb = IT'S — i^B,c and the conjugate ordering field m = With these variables the 

equation of state near the critical point is written as follows jsl: 

sfH = MV(B'"'^) ;t = ^ (1) 

In eq. (1) the scaling function f{x) is universal and properly normalized: /(O) = 1 and 
/(— 1) = Sj. Moreover, the function f{x) is regular at the point x = and follows a power 
law, f{x) ~ x^, for large x The constants B in eq. (1) are nonuniversal amplitudes 
and the critical exponents /5, 7, 5 are fixed by the 3d Ising-QCD universality class (/3 ~ |, 
7~|, 5~5). 

In the notation of reference [S] the critical equation of state for strongly interacting 

matter, given by the generic equation (1), is specified as follows [5|: 

B^"^t Bi\m\ . 

y = fix) ; X = , y = (2a) 

00 

/(x) = l + (x^O) (2b) 

n=l 



oo 



/(x)=x^^/rx-^"^ (x^oo) (2c) 

n=0 

For 3d Ising systems, the coefficients {fn \ f^) have been evaluated in reference [si] for 
n < 5, using renormahzation group techniques. The truncation of the series (2b, 2c) at the 
maximal order n = 5, leads to a good approximation of f{x) in the limiting regimes {x ~ 0, 
X oo) but also it guarantees a smooth transition from one region to the other around the 
point X ~ 1. 

Within the framework of eqs. (2) one may approach the critical point following distinct 
paths in the QCD phase diagram: for t = one moves along the critical isotherm, \pb\ = 
-Bcl^l^^*^, and the baryon-number susceptibility, xb = f^, develops a power-law singularity 
for m = 0, X-B ~ I'^l"^) a clear signal of a critical point (e = ^^). For t < 0, m = 
we move along the coexistence line, \pb\ = B\t\^, which is a limiting curve of the family 
(2a) when m ^ 0. In the general case (m ^ 0) eq. (2a) penetrates the crossover regime 
for t > with a smooth dependence of ps on t. Finally, in the limit a; ^ oo the system 
follows the critical isochor {ps ^0) along which the baryon-number susceptibility obeys 
the power-law xb = C'xKI with the amplitude = joo^-y/p ■ This landscape of universal 
power laws, extracted from exact equation of state, may be confronted with measurements 
of the observables: ub, ps, T in relativistic nuclear collisions, in order to capture the critical 
region of strongly interacting matter. 

For the net-baryon density hb we have used two different averaging scales in rapidity in 
order to verify that the nature of critical phenomena incorporated in the equation of state 
(2) does not depend on this scale. The ffist choice corresponds to the total rapidity length 
A available in the colhsion and leads to the representation: ^^(A) = AyjA'^J^{2sii\\i^)~^ 
where Ay^ is the number of wounded nucleons (participants). The other, rather extreme, 
choice is a scale of homogeneity at midrapidity {5y ~ 2) leading to the expression: ^^(0) = 
^niin i^Jl^ where j denotes the baryon number per unit of rapidity, in the central 
region. 

In Table [T] we summarize the experiments with nuclei (A+A') and the corresponding 
measurements of the relevant observables, in a wide range of energy and size of the colliding 
systems. The baryonic densities have been calculated using either the formula 6| 

B- B = (2.07 ± 0.05) (j9-p) + (1.6±0.1)(A- A) (3) 
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TABLE I: Measurements of the relevant observables in experiments with nuclei. 



A+A' 








/Xij(MeV) 


T[MeV) 


Ref. 


(S+S)cen. 


19.4 


54±3 


8.6±1° 


220±22 


180.5±10.9 


[8j,[9],[10j 


(S+Ag)cen. 


16.3 


90±10 


15±1'' 


242±18 


178.9±8.1 


MM 


(S+Au)cen. 


13.5 


113 


22±2'' 


175±5 


165±5 


[8j,[9],[l^ 


(Pb+Pb)cen. 


8.77 


349±5 


108.0±2.6^ 106.2±3.2'' 


381.6±6.7 


144.6±2.3 


[12j,[13j,[]J 


(Pb+Pb)cen. 


12.3 


349±5 


84.8±2.3^ 77.2±2.8'' 


296.0±6.4 


151.7±2.9 


[12j,[13j,[]J 


(Pb+Pb)een. 


17.3 


352±12 


67.7±7.6'' 


247.4±5.7 


156.1±1.6 


[6],[]J 


(Pb+Pb)o_5% 


17.3 


362±12 


82.7^ 


248.9±8.2 


157.5±2.2 


[7],[15j,[16j 


(Pb+Pb)5_i4% 


17.3 


304±16 


67.4* 


235.2±8.5 


150.6±3.2 


[7],[15J,[16J 


(Pb+Pb) 14-23% 


17.3 


241±16 


50.8'' 


223.7±9.9 


156.6±3.7 


[7],[15j,[16j 


(Pb+Pb)23-31% 


17.3 


188±16 


35.7** 


210±11 


154.7±4.2 


[a,[15j,[16j 


(Pb+Pb)3l_48% 


17.3 


130±14 


21.6* 


213±16 


153.2±5.9 


[a,[15j,[16j 


(Si+Si)cen. 


17.3 


se.i'^ 




253±11 


162.7±4.1 


M 


(C+C)cen. 


17.3 


24'= 




256±13 


166.1±4.3 


[14] 


(Au+Au)cen.— per. 


130 


381-4.9 


(0.090-0.018)"^ 


28.6-17.6 


165 


[ia,[18j,[19],[20 



"Calculated using A — A. 

'Calculated using — K^. 

'^It corresponds to the assumption A^j — 2A. 



or the formula 



a 



1 + 2b 



a = 1.07, b = 0.1 (4) 

In Fig. [1] the freeze-out states in experiments at the SPS and RHIC are located in the 
diagram nB(A) versus fiB- The experimental points belong to three distinct classes of 
measurements: (a) central and peripheral Pb+Pb collisions at the SPS (NA49 experiments), 
(b) central collisions of light nuclei (S, Si, C) at the SPS (NA35 and NA49 experiments) 
and (c) central and peripheral Au+Au collisions at RHIC. In the states of class (a) and in 
particular for the highest SPS energy {^/s ~ 17 GeV), the freeze-out temperature remains 
practically constant (Table [T]) and the actual value measured in these experiments is fixed, on 
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50 100 150 200 250 300 350 400 



Hb (MeV) 

FIG. 1: The baryon number density n^CA) versus baryochemical potential /i_B is illustrated for a 
series of freeze-out states. The solid line represents the best fit solution in the description of the 
Pb+Pb {^/s =17.3 GeV) central and peripheral data with the critical isotherm equation of state. 

the average, at the level (T) ~ 155 MeV. As a result we may assume that the experimental 
points of class (a) belong to an isotherm which, as shown in Fig. 1, develops a very large 
derivative |^ near the value — 200 MeV of the chemical potential. This observation 
implies that the baryon-number susceptibility of strongly interacting matter may become 
infinitely large in this region of the phase diagram (T ^ 155 MeV, — 200 MeV). It is 
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pq 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



■ Calculted using A 

O© Calculted using K 



8.77 GeV 



O 



Pb+Pb 
17.3 GeV 




23-31% 



nB.c(0) = 



1.53 Critical 
Point 



12.3 GeV 



(Pb+Pb) 



cen. 
i 17.3 GeV 



, ~S+Ag 



31-48% 




S+S 



Au+Au 
130 GeV 



I I I I I I I I I I I I I I I I 11*1 I I I I I I I I I I I I I I I I 



50 100 150 200 250 300 350 400 

Hb (MeV) 

FIG. 2: The same illustration as in Fig. [H using for 71^(0) the averaging scale at midrapidity. The 
solid line is produced through a best fit on the Pb+Pb {^/s =17.3 GeV) central and peripheral 
data (open circles). 

therefore suggestive that the freeze-out states in class (a) form a critical isotherm associated 
with the 3d Ising-QCD critical point. 

In order to verify quantitatively this hypothesis, starting from the above observations, one 
must employ the critical equation of state discussed in the first part of this Letter (eqs. (1) 
and (2)), in the case t = (critical isotherm). In other words, one must be able to describe 
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0.001 



150 160 170 180 190 200 210 220 230 240 250 260 270 280 

Hb (MeV) 

FIG. 3: The baryon-number susceptibility xb, for T = Tc versus hb- The soHd hue corresponds 
to the derivative of the best fit solution shown in Fig. [TJ The data points are produced using 
neighbouring (cycles) or next to neighbouring (squares) differences. 

the freeze-out states of class (a) in the diagram of Fig. [H with the equation of state: 



Ins (A) - nB,c\ = Be 



(5) 



With the normalization condition ns = for fiB = 0, the best fit solution [x^/dof ^ 0.16) 
is illustrated in Fig. 1 and corresponds to the actual values of the parameters: ub^c = 0.36, 
I^B,c-'2U MeV (5,~nB,c). 




140 145 150 155 160 165 170 175 180 

T (MeV) 



FIG. 4: The equation of state \pb\ versus T, for m = const. The curve m = 0, T < Tc corresponds 
to the coexistence Hne and for T > Tc to the critical isochor. The experimental freeze-out states 
are located within the area spanned by the family of curves of constant m. 

In Fig. [2] the same procedure is applied for a different averaging scale in rapidity, cor- 
responding to a baryon- number density 71^(0), as explained in the first part of this Letter. 
The critical solution is not affected by this change {^b,c — 218 MeV, /dof ^ 0.51) and 
the overall outcome of this treatment seems to be robust against the change of the averag- 
ing scale in rapidity. In what follows we use, therefore, ns(A) as a typical baryon-number 
density, on the average. 
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In summary, we have been able to isolate a class of freeze-out states corresponding to 
central and peripheral Pb+Pb collisions at ^/s ~ 17 GeV which form a critical isotherm 
associated with a QCD critical point at a location in the phase diagram fixed by the values: 
Tc ~ 155 MeV, fiB,c — 214 MeV. The baryon-number susceptibility, xb ~ \f^B — Ai_B,c|~*~, 
diverges at the critical point according to a universal power law, fixed by the critical exponent 
5 (eq. (5)). In Fig.[3l the singularity of xs corresponding to solution (5) is illustrated together 
with the results of a crude differentiation of nearby experimental points in Fig. 1. The errors 
are large but the trend is consistent with the derivative of the solution (5). 

In this last part of our treatment we attempt to exploit, beyond the critical isotherm, 
the universal properties of the critical equation of state which are rigorously incorporated 
in the representation (2). To this end, in Fig. H] we consider the diagram \pb\ versus T in 
which the freeze-out states are properly located together with the critical point found in our 
approach. With the help of eqs. (2) we have constructed, in the same figure, the family of 
curves corresponding to different values of m and representing the critical equation of state 
in its full capacity (critical surface). Across this family we have drawn the curve x — 1 
(eq. 2a) which, as it was already explained, corresponds to a smooth transition from the 
representation (2b) to the expansion (2c) of the equation of state. For m = 0, T < 
the coexistence curve has a normal scaling behaviour near the critical point, \pb\ ~ l^l'', 
and for B = 1 a. good description of the peripheral freeze-out states in class (a) is obtained 
along this curve. In general the lines of constant m, in Fig. HI provide us with a guidance 
for the identification of the experimental freeze-out states as critical states. In fact the 
cluster of experimental points in the area \m\ < 0.16, corresponding to a domain in chemical 
potential 214 < fis ^ 250 MeV, is an example of freeze-out systems consistent with the 
critical equation of state. In particular, the collision Si+Si {^/s ~ 17 GeV) freezes out 
close to the critical point (Fig. H]) and as a consequence, power-law density fluctuations are 
expected to occur in this process [2l|]. On the other hand, the freeze-out system at RHIC 
(Au+Au, ^/s ~ 130 GeV) violates strongly the critical equation of state and therefore a 
drastic change of the collision parameters is needed (lower energy, smaller systems) in order 
to capture freeze-out states consistent with criticality [22]. 

In conclusion, we have shown that the universal form of the critical equation of state for 
systems belonging to the 3d Ising universality class [5] can be fully exploited in experiments 
with nuclei in order to identify critical QCD states at the freeze-out regime. In this approach 
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we have found that existing measurements of ub, T, hb at SPS energies, reveal, within 
experimental errors, a class of freeze-out states belonging to the critical isotherm of the QCD 
universality class and leading, as a consequence, to a divergent baryon-number susceptibility. 
This phenomenon is a strong indication for the existence of a QCD critical point in the 
neighbourhood of these states {^b,c — 214 MeV, ~ 155 MeV). The overall analysis 
suggests that new, high precision measurements in experiments with light nuclei (S+S, 
Si+Si, C+C) are needed (SPS-NA61, RHIC at low ener gies ) in order to discover the exact 



location of the QCD critical point in the phase diagram 



2^,|23|. 
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